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ABSTRACT 


A theoretical solution of the equations of motion is developed for 
damped rigid rotors with harmonic forces at the supports. Support de- 
flection amplification and the phase angles between the displacements 
and the harmonic forces are computed. The harmonic forces represent the 
effects at the supports of unbalance forces whose planes of unbalance 
are assumed to have an equal distribution for the probability of their 
axial locations. Also, support stiffness and damping characteristics 
are assumed to be symmetrical. Under these conditions, which are repre- 
sentative of cylindrical roller bearings in rigid pedestals, a low-speed 
resonance appears when the center of gravity does not lie on the mid- 
plane. Other calculations show that even with a midplane center-of- 
gravity location, a low-speed resonance occurs for slender rotors but 
not for disks. Deflection amplification at low-speed resonances is 
generally smaller than at the major critical speed. 
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DEFLECTION ANALYSIS OF DAMPED ROTORS WITH HARMONIC SUPPORT FORCES 


by Richard • H . Cavicchi 
Lewis Research Center 


SUMMARY 

A theoretical steady-state analysis was performed to study damped 
vibrations of flexibly-mounted rigid rotors subjected to harmonic forces 
at the supports. The rotor-support system is mounted on a firm founda- 
tion. Viscous damping is assumed in which the viscous force is propor- 
tional to the rate of support displacement. Also, the support stiffness 
and damping characteristics are assumed to be symmetrical; symmetrical 
stiffness is a valid approximation for the characteristics of cylindrical 
roller bearings mounted on rigid pedestals. The harmonic forces are rep- 
resentative of the effects at the supports of unbalance forces. Thus, 
there is an equal distribution of the probability of the axial location 
of the planes of unbalance. A proportional relation is assumed between 
the force amplitudes at the two supports. The results of this analysis 
are calculated values for the support deflection amplifications, and the 
phase angles between the displacement and the harmonic driving forces . 

Seven illustrative cases are selected to display the results of 
the analysis graphically over a wide range of rotational. speed. The 
parameters varied in these cases are damping, rotor shape (which is 
varied from a thin shape to a disk) , center-of -gravity location, and 
force ratio at the supports. 


The graphical results reveal that in addition to resonance at the 
major critical speed, a lower-frequency resonance occurs when the center 
of gravity is not located at the midplane. Other calculations show that 
even. with the center. of gravity at the midplane, a lower-frequency reso- 
nance occurs if the polar-to-diametral moment-of-inertia ratio is less 
than 0.7. Thus, disks should not experience the low-frequency resonance. 
Deflection amplification is generally greater at the major critical speed 
resonance than at the lower-frequency resonance; however, the low- 
frequency amplification may be. greater for highly damped rotors . 

Results from varying the center-of-gravity location suggest that 
there are Certain locations that should be avoided in design because of 
potentially large deflection amplification. 


INTRODUCTION 

The problem ofi vibration effects in rotating machinery has been in- 
vestigated for more than a century. Understanding of this problem has . 
advanced considerably in, recent years as the need for high-speed rotors 
has intensified the effort in rotor dynamics . Numerous investigators 
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have used various models, approaches, arid conditions in this work. Among 
various aspects of this problem that have been studied are critical 
speeds, deflections, dynamic stability, transient effects, and the ef- 
fects of acceleration. 

The ability to predict and control rotor dynamics , as reflected 
by these aspects, is vitally important in the design of modern high-speed 
rotating machinery. Such machinery is frequently designed to operate at 
speeds above the rigid-body critical speeds. Thus, in startup and shut- 
down, rotors must pass safely through these critical-speed regimes. 
Usually, acceleration or deceleration is performed as rapidly as possible 
in order to avoid a damaging buildup of resonant amplitude.. As improved 
analytical and design techniques are developed to predict and control 
critical speeds and resonant amplitudes, 5 the designer will have greater 
freedom to design for operation near' or even at critical speeds. 

Reference 1 presents an analysis that identifies critical speeds 
for vibrations of undamped rotors on firm foundations . The critical 
speeds include major critical speed (precession frequency equal to rotor 
frequency) and nonsynchronous critical speeds (other speeds at which de- 
flection amplification may occur). Yamamoto’s analytical and experi- 
mental work, reported in reference 2, shows that bearing defects cause 
nonsynchronous critical speeds. In. ref erence 3, Gunter inquires into 
the dynamic stability of rotor-bearing systems. 

The analysis in reference 4 locates critical speeds for free vibra- 
tions of undamped rotors on flexible foundations. It studies the effects 
of variations in rotor speed, rotor and foundation shapes , and rotor-to- 
foundation spring-constant ratio. A comparison of the firm- and 
flexible-foundation results of references 1 and 4 appears in reference 5. 

In reference 6, Gunter and DeChoudhury provide ALGOL computer pro- 
grams to study steady-state effects of rotor unbalance, dynamic sta- 
bility, and transient response . The linearized equations used in this 
reference contain cross-coupling stiffness and damping coefficients. 

The present report extends the analytical work on rotors of refer- 
ence 1 to include viscous damping and harmonic disturbing forces. It 
determines critical speeds, support deflections , and phase, angles in re- 
sponse to harmonic equal and unequal forces at the supports. The har- 
monic forces at the supports represent the effects at the supports of 
unbalance forces that have ari equal probability for axial location or 
magnitude. Another application of the harmbnic forces in this report is 
to the magnetic forces in a turboalternator. Effects of this type of 
force on a rotor are studied in reference 7. Experimental values of 
cross-coupling coefficients are usually difficult to determine. The 
present report therefore uses a simplified version of the equations of 
motion that omits cross-coupling terms. The simplified analysis is used 
in this report to show trends resulting from varying damping, force ratio 
at supports, center-of -gravity? location, and momerit-of -inertia ratios. 
Seven combinations of parameters are- Used to “illustrate the trends. 
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DESCRIPTION OF ROTOR-SUPPORT MODEL 

Figure 1 is a sketch representing the support system of a rotor 
with damping. This study considers lateral motion (in the y- and z- 
directions) but no axial motion (x-direction) . The symbols are defined 
in appendix A. 

The support spring constant k/2 is assumed to have the same value 
for both supports in both the y- and z-directions . This is a good 
approximation for cylindrical roller bearings mounted in rigid pedestals. 
Likewise, a common value of damping factor b/2 is used at both supports 
in both directions. Static deflections are neglected; this would apply 
to conditions in space. 


Harmonic Exciting Forces 

Forced vibration due to the unbalance of eccentricity can be repre- 
sented by a harmonic, force acting at the center of gravity and a harmonic 
couple. The amplitudes of this force and couple are independent of time. 
These amplitudes are, however, functions of rotational speed which is in- 
dependent of time in the steady-state solutions of this analysis. 


For convenience in this study, the resultant force at the center of 
gravity is resolved into lateral components acting at the supports. The 
couple is replaced by the equivalent effects, of the moments about the 
center of gravity of the force components at the supports. 

The left support in figure 1 is designated as support 1 for identi- 
fication and the right support is support 2. The force amplitude at sup- 
port 2 is assumed to be proportional to the value at support 1. In the 
examples given in the report the proportionality factor r has been 
assigned values of 0, 1, and 2 to represent possible exciting force con- 
ditions due to unbalance. 

Viscous damping at the supports is assumed in which the retarding 
force is taken to be proportional to the rate of displacement at the 
shaft supports. 


Rotor Shape 

. Rotor configurations that vary from shapes as thin as pencils to 
the other extreme represented by disks are considered in this analysis. 
Shape is expressed in terms of the rotor polar moment of inertia Ip and 
its diametral moment of inertia I D . The analysis is made nondimens ional 
for generality; therefore, these two moments of inertia are expressed as 

tt 2 = I d /M£ 2 


( 1 ) 
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and 


* 3 - y^ <2) 

Both parameters have small values for slender shapes. For example, 

TT 3 = 0.0149 for a rod whose length is 10 times its diameter. Its value 
for a thin disk is 2. The parameter called the disk effect, is 

zero for a concentrated mass and is infinite for a disk having all its 
mass distributed over a large radius. 


Center-of-Gravity Location 

Although the center of gravity is assumed to remain on the shaft 
centerline, its axial location is varied. In the examples, four typical 
values are specified for this location, from midplane to outboard quarter 
point, the last of which represents an overhung configuration. 


ANALYSIS 
Exciting Forces 

The disturbing force components. F y and F z in the model of 


?1 

figure 1 are assumed to be harmonic and represented by 

F = F (to) cos tot 


(3) 


and 


F = F (to) sin cot (4) 

acting at support number 1. The corresponding components at support 2 
are assumed to be the multiples 


sb r V 

y 2 771 


and 


F 


z 2 


r z F Zl 


(5) 


( 6 ) 


The force amplitude F(to) is assumed to be independent of time as dis- 
cussed earlier. In particular, F(co) — moo^e represents an unbalance 
force due to eccentricity e. On the other hand F(oo) = a constant rep- 
resents the magnetic forces in a turboalternator as shown in reference 7. 
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Moments About Center of Gravity 

The resulting moment about the y-axis through the center of gravity 
is 


’TCy " -»l r z 1 + V) 

The moment about the z-axis is 

= *1^ - *2 F y 2 (8) 

When equations (3) to (6) are used, the two moment equations become 

Ttly = .F(A 2 r z ~ ^ sin wt (9) 

= -FCfl^y - £i)cos cot ( 10 ) 


Equations of Motion 

The following equations of motion are a mathematical approximation 
of the rotor-support model of this study. 


*2^1 , W , b .. , . , , k 
~~i~ + ~ F~ + 1 < y i + y 2> + 2 'n 


+ + J (^ + y 2 ) + 2 ^1 + 3^ " + r y) cos ut ( 13L ) 

£ Mz & Mz 

■ --£ + — J - 2 - + \ (z 1 + z 2 ) + j (z 1 + z 2 ) = F(co) (1 + r z )sin cot (12) 


p / • .. 13 / .. *.\ b • , b „ . k rt . k . 

% y 2 " y l^ ” a Z 2 " Z l ) - 2 l 2 z 2 + 2 A l z l " 2 £ 2 Z 2 + 2 £ 1 Z 1 


F(co)(£ 2 r z - £^)sin cot 


(13) 


p * • • \ i -0 /,, aa . . . b . . b . . . k . ft k . 

I (z 2 - z x ) + % (y 2 - y x ) + 2 £ 2 y 2 ■ 2 £ l y l + 2 £ 2 y 2 “ 2 £ l y l 


= -F(co)(A 2 r^ - £^)cos cot (14) 


Laplace Transforms 

It is convenient to take the Laplace transforms of the equations of 
motion. The following development assumes that the support displace- 
ments and derivatives at time zero are zero. 


f 
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With the following type of notation for a Laplace transform 

[y(t) ] = Y(s) - Y (15) 

and zero initial conditions , the equations of motion can be expressed as 

mi, 

L i\ i 


& 


k’ 


2 X. b . 

s + Y s + 2 


Mil 


2 2 , b 


1\ l 


mi 


1 2 , b 


-Y 

1\£ 


s + Y 


I D 2 M 1 kZ l\ 

r 3 + — s + T - 


2 \& 


! ) + h 


s(r y + l)F(w) 

(16) 

2.2 
s + U) 

u(r + l)F(rn) 
z 

(17) 

2 , 2 
S +0) 

, f 1 D 2 b£ 2 

'2\£ S + 2 

k£ 2 ’ 
S + 2 


m(£ 0 r - ) F (co) 

Z Z 1 


2 2 
s + w 


(18) 




I. biL kJ l, 

D 2 1 . 1 

— s + — r— s H — 


1-is O btp k£_ 

+ hlr s + -r B *^'- z 


) - •) * h(? •) 


-s(^2 r v “ £ x) F (“) 


2 , 2 
s 4* to 


(19) 


The coefficients of the Laplace variables are constants A-^, A 2 , B-^, B 2 , 
and C as defined in the symbols list. Use of these constants in equa- 
tions (16) to (19) yields 


A 2 Y 1 + A 1 Y 2 


s(r y + 1 )F(uj) 

2,2 

s + to 


( 20 ) 


A 2 Z 1 + A 1 Z 2 = 


oi(r + l)F(co) 
z 

. - 

S + ft) 


( 21 ) 


~ CY ! + CY 2 + B 1 Z 1 - B 2 Z 2 


w( V z " V F(u>) 

2,2 

S +0) 


( 22 ) 


-s(£ r - £ )F(u) 
-B,Y. + B Y, - CZ, + CZ„ = s- 2 


11 2 2 


2,2 

S +0) 


(23) 
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Deflection Solutions 

The object of this analysis is to- solve equations (20) to (23) for 
the deflections represented by and Y 2 at each support. Due to 

symmetry, the deflections Z]_ and Z£ differ from Y]_ and Y£ only 
by 90° in phase angles; therefore, no solutions will be developed for Z^ 
and Z£. The solutions for Yj_ and Y 2 are given by 


and 


where 


Y l“ 


F(to) 

2 , 2 
'vS + oj , 


N, 


(24) 


Y 2 = 


F(m) ' 

2 2 
+ 0 ) j 


N, 


(25) 


N l = 


N 2 " 


(r 

y 

+ l)s 

0 

0 


+ 1)0) P 

A 2 

A 1 

u 2 

r - &,)o) C 

Z X 

B 1 

- B 2 

l~^ £ 2 r y £ 1^ S B 2 

-C 

c 

A 2 

(r + l)s 

y 

0 

0 

0 

(r z + l)w 

A 2 

A 1 

-C 

( Vz " £ l )u 

B 1 

-B 2 

1 — 1 

PQ 

1 

~ ( ' SL 2 V y “ V S 

-C 

C 


a 2 a 1 0 

0 


) « 

0 0 a 2 

A 1 



-C C B x 

- B 2 



" B 1 B 2 - C 

c 



(26) 


(27) 


(28,) 


Expressions for N-|_, N 2 , and D obtained from expanding the deter- 
minants in equations (26) to (28) are given in appendix B. 
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All variables are made dimensionless in this analysis and are sum- 
marized here for reference. 


s 

^ /k/M 



( 29 ) 


( 1 ) 

( 2 ) 


TT, 


7T 


y k/M 

_ _L_ 
10 " M 

L -*1 
L 1 A 


J 2 a 


( 30 ) 

( 31 ) 

( 32 ) 

( 33 ) 


Deflection Amplitude and Phase Angle 

This report studies the frequency response to a sinusoidal input. 
Mathematically, this is achieved by setting s = i . In dimensionless 
form, this is 


n, = in. 
1 4 


( 34 ) 


Then 


Y _ a + ic 
F(m) b + ih 
k 


from which the magnitude is 

fj 77 

Ms 2 ** 2 



( 35 ) 


and the phase angle between the deflection and exciting force is 


( 36 ) 
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0 


tan 


-1/ eg - ah \ 
\ag + eh/ 


(37) 


Detailed expressions for a, c, g, and h are presented in appendix B. 


Normalized Deflection 

The ratio in equation (36) is normalized by dividing by its value 
at zero rotational speed (when ir^ = 0) which represents the deflection 
at zero exciting frequency or large stiffness; therefore, equation (36) 
is expressed by the amplification factor 


V 



Y 1 



F(u)/k 


Y 1 



F(io) /k 

^4 

=0 


, 2 . 2 
i a + c 

2 j. v2 



(38) 


The normalizing factors are derived in appendix B. Then, by equa- 
tions (B16) and (B17) , equation (36) becomes 



Y 1 


F(oo)/k 

Y 1 

o 

II 

F(oj) /k 



2[L 2 (2r + 1) - L x ] 


(39) 


for support 1. The corresponding expression for support 2 is 


y 2 = 



Y 2 


F(ai)/k 

Y 2 

o 

II 

F (to) /k 


, 2 , 2 

a 2r + c 2r 

2 , ,2 
_g_ + h 


2[r(L 2 - L 1 ) - 21^] 


(40) 


CATALOG OF EESULTS 
Parameter Selections 

An abundance of parameter combinations is possible for making cal- 
culations in this analysis. Seven particular combinations were chosen 
to display the effects of parameter variations. The intention was to 
provide meaningful results with minimum data production. 
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The seven parameter combinations selected are tabulated in table I. 


Data Plots 

For each of the seven cases of table I, calculations were made for 
the amplification factor y (eqs. (39) and (40)) and phase angle 0 
(eq. (37)) for both supports. Data plots were made in which y and 0 
are each plotted against rotor rotational parameter tt^. Each point on 
the data plots represents a steady-state point. 

A negative angle represents a phase lag between the driving force 
and the support displacement. 

Case 1 . - The case 1 data plots in figure 2 are typical of the fre- 
quency response of classic second-order systems presented in textbooks. 
Perfect symmetry is assumed for case 1. That is, the center of gravity 
is at the midplane (L]_ = L 2 = 0.5) and both supports are subjected to 
forces of the same amplitude (r_ = r z = r ■ 1) . Figure 2(a) shows the 
effect of damping parameter on amplification factor and figure 2(b), 

on phase angle. 

The amplification factor in figure 2(a) is very sensitive to the 
amount of damping in the region about ir^ « 1. This, of course, is the 
condition for which the rotational frequency equals the natural fre- 
quency. Because the amplification factor becomes infinite in an. undamped 
system at ir^ = 1, it is also referred to as the major critical speed. 
Farther removed from the condition of ir^ = 1, the effect of damping on 
amplification factor is minor. 

Figure 2(a) also shows the typical effect that the resonance peaks 
occur at gradually lower rotational speeds as damping is increased. No 
peak at all is reached for values of tt^q greater than a value of 
about 1.2. 

Forced vibrations with viscous damping are harmonic and have the 
same frequency as the driving force. The lateral support deflections, 
however, are out of phase with the driving force. The phase lag is pre- 
sented in figure 2(b) as a function of damping and rotational speed 
parameter. At the major critical speed, all systems show a phase lag of 
90° irrespective of the amount of damping. Below this value, highly 
damped systems have the greatest phase lags. Above = 1, the opposite 
is true. At rotational speeds near the major critical, the phase angle 
is extremely sensitive to rotational speed for lightly damped rotors. 

Figure 2 shows that y approaches 1 and 0 approaches zero for all 
values of damping as the rotational speed approaches zero. Thus, at low 
rotational speeds the deflection is the same as though the force were 
applied statically. Also, at high rotational speeds the phase lag for 
all values of damping parameter approaches 180°. 
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Case 2 . - The plots in figure 3 show the effect of damping for the 
extreme condition in which support 2 experiences no exciting force 
(r = 0) . The ordinate in figure 3(a) is given as a deflection-speed 


parameter 



A normalizing factor at zero rotational 


speed is not used for support 1 for case 2 because it is zero for r = 0 
and = L 2 = 0.5 (see eq. (B16)). All other deflection parameter plots 

in this report are normalized except for figure 3(a). Figure 3(b) shows 
the amplification factor for support 2. Solutions for equations (20) to 
(23) for Zj_ and Z 2 were obtained, although not presented in this re- 
port. They reveal zero deflection in the z-direction at support 1 at 
zero speed when r is zero. 


Equations (B18) and (B19) relate r and the center-of -gravity lo- 
cation for zero rotational speed and zero deflection at each support. 

When the center of gravity is at the midplane (Lj_ = L 2 ) , equation (B18) 
shows that r = 0 for zero deflection at support 1. But, according to 
equation (B19) there is no center of gravity location for which r = 0 
when the deflection at support 2 is zero at zero rotational speed. 

Figure 3(b) shows that there is no deflection amplification for sup- 
port 2 when the damping parameter tt^q is above about 0.6. By compari- 
son, the critical damping parameter is 2.0 for free vibrations with mid- 
plane center of gravity. Critical damping may be defined as the value 
for which maximum amplification factor is 1. 

In the phase angle plots of figures 3(c) and (d) all curves do not 
pass through 90° at = 1 as occurs for case 1. This results from 
the lack of symmetry of the, exciting force. 


Case 3 . — In. case 3, support 2 is subjected to twice the force of 
support 1 (r = 2). Figure 4 presents the curves for this case. Fig- 
ure 4(a) shows a critical value of tt^q at a little greater than 0.8 
for support 1. The corresponding value for support 2 is 2.2 (fig. 4(b)). 

In the phase-angle plots of figures 4(c) and (d) , the tt^ = 1 
crossings are more closely bunched than for case 2. Otherwise, the 
curves in figures 4(c) and (d) are similar to those in figures 3(c) 
and (d) . 


Case 4 . - Figure 5 presents plots that show the effect of center- 
of -gravity location. When the rotor center of gravity is located away 
from the midplane, an additional set of resonances appears. This set of 
resonances is also influenced by the support stiffness and damping char- 
acteristics, which are assumed to be symmetrical in this analysis. The 
result is that with a distribution of equal probability of the axial lo- 
cation of the planes .of unbalance, .the probability of exciting a low- 
frequency resonance increases if the center of gravity is not located 
at the midplane. The rotor rigid body cylindrical and conical modes 
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are coupled together when the center of gravity does not lie on the mid- 
plane . 

Figures 5(a) and (b) show a resonance peak occurring at a rotational 
speed below the major critical for each curve except the = 0.5 curve. 
The low frequency resonances are much smaller in both amplification and 
magnitude than those at the major critical speed. This means that a 
given amount of damping is more effective in displacement attenuation at 
the low-speed resonance. 

Figures 5(a)' and (b) show that as the center of gravity is located 
progressively farther from the midplane, the resonance peaks occur at 
lower rotational speeds. This observation is true for both the low- and 
high-frequency sets. 

A conclusion from figures 5(a) and (b) is that even though addi- 
tional resonances may occur when the center of gravity location is 
varied, the more severe deflection amplification occurs at the major 
critical speed. 

Comparison of figures 5(a) and (b) reveals that at support 1, maxi-: 
mum deflection amplification occurs when the center of gravity lies at 
the quarter point (L^ = 0.75). At support 2 the worst condition occurs 
when the center of gravity is at the midplane. Furthermore, greater 
amplifications occur at support 1 than at support 2. 

Case 5 . - In case 5 the center of gravity is located at the support 
subjected to the larger driving force (support 2). The damping param- 
eter varies for this case which is presented in figure 6. The curves 
for = 0.1 are repeated from figure 5. The two sets of resonance 

peaks are well defined. At both supports, the highest resonance peak . 
of the low-frequency set occurs with a damping parameter of 0.5. Fig- 
ure 6(a) shows no low-frequency resonance peak at zero damping. Instead, 
there is an antiresonance or minimum amplification. At support 2 with 
zero damping there is no resonance or antiresonance among the low- 
frequency set. On the other hand, the low frequency peaks in figure 6(b) 
are higher for tt^q >. J 2 than the high frequency peaks in figure 6(a). 

The critical value of tt^q for support 1 is about three for the 
low-frequency set as shown in figure 6(a) and 1.8 for the high-frequency, 
set. At support 2, figure 6(b) , the corresponding values for tt^o are 
2 . 8 and 0.6. 

Comparison of figures 6(a) and (b) suggests that deflection ampli- 
fication is more severe at support 1 than at support 2 when the center of 
gravity and the larger force are at support 2. 

Case 6 . - The results for case 6 appear in figure 7. This case pre- 
sents the effect of the polar-to-diametral moment-of-inertia ratio TTg 
on deflection amplification. At both supports the high frequency peaks 
all occur at tt^ = 1. In figures 7(a) and (b) , no low-frequency peaks 
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occur for values of ^3 greater than 1. Therefore, disks (113 = 2) ex- 
perience large deflection amplification at only the major critical speed 
whereas slender shapes have two resonant speeds. The value of ^ 3=0 
is fictitious and merely sets a lower boundary to the calculations. The 
curves representing H 3 = 1 depart from the usual pattern at high rota- 
tional speeds. 

The phase angle plots in figures 7(c) and (d) show all curves pass- 
ing through 90° at = 1. This results because of the low value of the 

damping parameter = 0.1). The curves would not all pass through the 

90° point at higher values of damping parameter as figures 4(c) and (d) 
imply . 

Comparison of figures 7(a) and (c) shows that the phase-lag peaks 
of figure 7(c) occur at rotational speeds between the low-frequency reso- 
nances and the antiresonances of figure 7(a). This is also true for sup- 
port 2 (figs. 7(b) and (d)). 

Case 7 . - Figure 8 presents the variation of the disk effect over a 
range from zero to 10. No additional resonances occur at either support. 
It should be realized that the value of ^3 used in figure 9 is 2. Fig- 
ure 7 has shown that no low-frequency resonances occur when ^3 is one 
or greater; therefore , curves drawn for case 7 with it 3 values of 0.5 
and lower may be expected to show low-frequency resonance peaks. But, 
the result of figure .8 is that varying the disk effect does not re- 

sult in an additional resonance peak. The = 0 curve has no physical 
significance. 

All the phase-angle curves for both bearings pass through 90° at 
the major critical speed in figures 8 (c) and (d) . 


INTERPRETATION OF RESULTS 

Cross plots were made of some of the amplification curves to aid. in 
interpreting effects of the parameters varied. All amplification factors 
discussed in this section are peak amplification factors. They were read 
at the appropriate peak. 


Damping 

Figure 9 presents peak amplification factor as a function of damping 
parameter. Curves are shown for both supports. Cases 1, 2, and 3 are 
represented in figure 9(a) and case 5 results are in figure 9(b). 

Center of gravity at midplane . - The center of gravity in fig- 
ure 9(a) is located at the midplane. In this figure it is immediately 
apparent that all curves have the same slope on the logarithmic plot. 

They can all be represented by an expression of the form 
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yir 10 = constant y > 1 

which is linear on a logarithmic plot. The condition y = 1 corresponds 
to critical damping. 

In figure 9(a), the line for r = 1 
line for r = 0 as long as y > 1. The 
drawn for Yp/F(co) /kj whereas the r = 0 
factor y. 

In a comparison of the dashed set of lines (r = 0) with the solid 
set (r = 2) in figure 9(a), it is evident that the amplification factor 
for both supports increases as the force ratio r is increased. 

When r = 0, support 1 experiences greater amplification than sup- 
port 2. The opposite is true when r = 2. 

Center of. gravity at support 2 . - The center of gravity is located 
at support 2 for the model of figure 9(b). Figure 5 has shown that an 
additional set of resonance peaks occurs when the center of gravity is 
displaced from the midplane location. The low-frequency set is dashed 
in figure 9(b) and the high-frequency set is solid. 

The high-frequency set shown in figure 9(b) has slopes essentially 
parallel to the lines in figure 9(a). 

Support 1 is subjected to greater amplification than support 2, in 
figure 9(b), for both resonance sets. In figure 9(a) for r = 2, the 
opposite was observed. 

Comparison of the dashed and solid curves in figure 9(b) shows that 
for highly damped rotors, the low-frequency amplification exceeds that 
of the high-frequency set. In all such cases, however, the amplification 
factor is less than 2. 

The low-frequency set of curves in figure 9(b) reaches a peak at a 
damping factor of about 0.4. This implies that the amplification de- 
creases for the low-frequency set as damping factor is decreased below 
0.4. This circumstance is of no practical advantage, however, because 
of the large increase in amplification of the high frequency set. 


coincides with the support 1 
r = 1 line, of course, is 

line represents the normalized 


Center-of-Gravity Location 

Figure 10 presents the effect of center-of-gravity location on peak 
amplification factor as determined for case 4. Two sets of resonances 
appear in this figure for each support. 

It is appropriate here to discuss the relative magnitudes of the 
curves for the two bearings at a given center-of-gravity location. Over 
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most of the range shown, the amplification factor for support 1 exceeds 
that of support 2 despite the fact that support 2 experience twice the 
force of support 1, This is not a paradox as can be seen by noting equa- 
tions (B16) and (B17) . For all values of Lp and L2 except when they 
are equal, the magnitude given by equation (B17) exceeds that of equa- 
tion (B16); therefore, the normalizing factor used in calculating y for 
support 2 exceeds that for support 1. Thus , although the actual deflec- 
tion at support 2 is greater than that at support 1, its amplification 
factor y is less. 

The most important information given by figure 10 is that there are 
certain center-of -gravity locations that should be avoided in order to 
prevent large deflection amplification. For the conditions of case 4, 
locations given by Lp in the vicinity of 0.83 should be avoided. Fig- 
ure 10 shows peak amplification factors at support 1 at this location 
for both the low- and high-frequency resonances . Support 2 peaks occur 
at Lp = 0.6 for both resonance sets so these locations should also be 
avoided. Calculations were made at very small increments in Lp in 
order to define the location of these peaks accurately. The location of 
the peaks would be different fox other parameter combinations, of course. 

The amplification factor for support 1 exceeds that for support 2 
over most of the center-of -gravity range studied. Below about 0.61 and 
above about 1.2 for Lp, the support 2 amplification is greater. The lo- 
cation at which the amplification factor of the high-frequency resonance 
set ‘is the same for both supports is very close to the peak for sup- 
port 2. 


Polar-to-Diametral Moment of Inertia Ratio 

Figure 11 presents the effect of polar-to-diametral moment-of- 
inertia ratio on amplification factor under the conditions of 

case 6. The most prominent aspect of figure 11 is the small effect that 
7T3 has on amplification factor of the high-frequency set. 

For the low-frequency resonance set, in general, amplification fac- 
tor decreases as it 3 increases. Above about 0.65 for ^3, no low- 
frequency set exists. 

Over most of the range in figure 11, the amplification factor at 
support 2 is twice that at support 1. 


Support Force Ratio 

The effect of support force ratio r on amplification factor 
appears in figure 12. Damping parameter values of 0.1, 0.2, 0.5, 0.8, 
and 1 are used in this figure.. Cases 1, 2, and 3 were used in produc- 
ing figure 13. 
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At support 2, represented by dashed curves, the amplification factor 
increases approximately linearly with r. The support 1 curves decrease 
with increasing values of r. At r = 1 both supports experience the 
same amplification due to symmetry. As r approaches zero, the support 
1 amplification approaches infinity, as verified by the fact that equa- 
tion (B16) yields zero. The zero magnitude of equation (B16) occurs only 
when L 2 , of course. At r = 0, the amplification factor at sup- 

port 2 remains finite (see eq. (B17)). 

The reason that the amplification factor for support 1 is greater 
than that for support 2 is explained by the relative magnitudes of equa- 
tions (B16) and (B17) . Although the amplification factor for support 1 
is approaching infinity as r approaches zero in figure 12, the deflec- 
tion is actually approaching zero. Figure 3(a) verifies this obser- 
vation. 


CONCLUSION 

This deflection amplification analysis of damped rotors shows that 
in addition to the major critical speed, a lower-frequency resonance 
appears if certain design conditions exist. These conditions include 
symmetrical support stiffness and damping characteristics and equal 
probability for the distribution of the axial location of the planes of 
unbalance along with either (1) center of gravity not located at mid-? 
plane or (2) ratio of polar-to-diametral moment-of-inertia ratio less 
than 0.7; therefore, the low-frequency resonance should appear for 
slender rotors but not for disks. 
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APPENDIX A 
SYMBOLS 


k (Vi + ^ 10 V 2 + 1/2 ) 
k ( L 2 Tr i + Wl /2 + 1 ^ 2 ) 


real part of numerator in eq. (35) for support 1; see eq. (B4) 

real part of numerator in eq. (35) for support 2; see eq. (B6) 

. ( 2 L i*io*i L i\ 

■ “ ^ 2*1 + + Tj 

- +— — *Tj 

linear damping coefficient, kg/sec 

U V 2 V4 


C 1 

imaginary 

part 

of 

numerator 

in 

eq. 

(35) 

for 

support 1, 

c 2 

imaginary 

part 

of 

numerator 

in 

eq. 

(35) 

for 

support 2; 

D 

denominator of 

deflection solution^ 

see 

eq. 

(22) 


eccentricity causing unbalance 

amplitude of harmonic driving force, nt 

real part of denominator in eq. (35) 

imaginary part of denominator in eq. (35) 

2 

rotor diametral moment of inertia, kg-m 

2 

rotor polar moment of inertia, kg-m 

/=! 

2 

linear spring constant, kg/sec 


1 £ 

l .fa 

L 2 Z 

5L Laplace operator 

Z distance between bearings, m 
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M rotor mass, kg 

°)7U moment, nt-m 

N numerator of deflection solution; see eq. (22) 

r bearing force ratio; see fig. 1 

s Laplace variable, sec ^ 

t time, sec 

x,y,z rectangular coordinates defined in fig. 1 

Y Y(s) = <£[y(t)] 

Z Z(s) =£[z(t)] 

0 phase angle, deg 

y deflection amplification factor 

s//k/M 

2 

tt 2 I d /M£ = disk effect 

TT 0 1 /I*-. = rotor polar-to-diametral moment-of-inertia ratio 

j p D 

it £ w/^k/M = rotational-speed parameter 

tt^q b/v/kM = damping parameter 

to rotor rotational frequency = harmonic driving force frequency, 

rad/ sec 

Subscripts : 

P peak 

r condition for r = r = r 

y z 

y y direction; see fig. 1 

z z direction; see fig. 1 

supports 1 and 2; see fig. 1 


1,2 
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Superscripts : 

• first derivative with respect to time 

• • second derivative with respect to time 
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APPENDIX B 
ANALYSIS DETAILS 
General 

Expansion of the determinants in. equations (24) to (26) yields 
— — 

N 1 " _(r y + 1)s c2(A 1 + V + B 2 (A 1 B 1 + A 2 B 2 ) " (r z + " V 

w 

+ A 1 [U 2 r z - Jl 1 )wC(A 1 + A 2 ) - a 2 r y - J^MA^ + A 2 B 2 ) ] (Bl) 

N 2 " ”< r y + Ds[c 2 (A 1 + A 2 ) + B 1 (A 1 B 1 + A 2 B 2 )j + (r^ + 1)UA 2 C(B 1 - B^ 

- A 2 [(Jl 2 r z - Jl 1 )coC(A 1 + A 2 ) - (A 2 r y - ^(A^ + A 2 B 2 ) ] (B2) 

D = -C 2 (A 1 + A 2 ) 2 - (A 1 B 1 + A 2 B 2 ) 2 (B3) 

The values of a and c for Y^ in equation (35) differ from the 
values for Y 2 

For support 1 (see fig. 1), 
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For support 2, 
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) sf 2 ^1 + ^2 ^l\ 3 V 2 ^l(^l ^ 

c 2 = _(r y + 1) V4r 7r 2 7r 3 + 2 + *2 + TF 2 " *4 F + 2 


L it ( 3 L 

+ — + v 2 +ir 4b 


l)\o + (r z + + A) 


( L 1 - L 2 )tt 10 


7 r 2 3 L 2 


/ 5 L 2 ( L 1 + L 2) . , . *2 3 L 2 ( L l + L 2 ) 

+ (r y L 2 L-l)^ 4 2 + L 2 ir 2 + 2 V' 2 + 2 

+ T 2, + ’2 + ” 4 (ti’l 0 + (r z L 2 - V Vs'lO^ + f) ’ ”4] < B7 > 


For each support, the denominator D has components 


ki - M - ”4 + "iJ 


2 2 

4 . Tt^Tt -f 
2 10 + 


^ 


( L 1 + L I )’ i 2 C 


2 | VlO 2 I M+'-lf T U 

1 + * 10 ) + 2 + 2 ”2 + *4 4 + 16 
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Restriction: r = r = r 

y z 

All the calculated results presented in this report were obtained 
for the condition 


r = r = r 

y z 


(BIO) 


Therefore, the values of a and c as modified by this condition are 
presented below. 


Support 1 . 


lr 



. h( L i + L z) 

“r > 


- Al + "i - *) + 


(”io + x ) 


(A + A) 


ITT 7T TT 

2 2 10 

2 4 


+ 1 - IL, + 


’ioV L i - L 2 


2 3 VI 2 


) 


, 3L 2 1T 10 V.2 
8 2 


f L l " L 2 '\ 

^2^3^ 4 / 


- IT. 


+ (rL 2 - + « + 4 


2 /V L 1 + 2 


, T t !io \ V10 

+ L 1 4 + n 2/ + 2 


( ir 3 + ^ + 1T 2 7r 3( L l + 2) + 2 


“ It, 


L 1 L 1 + L 2 3 ^10 *2 > 1 - . . 1 

T + 4 + 1~~ + T (ir 3 + 1) + 8 


(Bll) 
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c lr = (r + 1)K 


-"zH - *) + - 


4 T JhJLti) 

2 + 2 L 1 3\ 2 ) 


\ ^2 
2 2 


( L 1 + L z ) 


+ vL + l 

+ 8 + L 2 


7T — 7F TT 7T l * 
2 2 3 10\ 


L 1 " V 


-.(4 


10 


r 

L , 


+ ( rL 2 - Ll )< 


( L J + L l) 


+ + — + * 2 


* 3 ( l i + l) 


2 2 2 

L 1 L 1 + L 2 *10 

— + 5 + ~ 8 ~ + 1 T 2 (ir 3 + 1} 


^(Drio 


(B12) 


Support 2 . - 


a 2r " " (r + US'lvl 


4 2 4 ill ( L 1 + L 2r2 ^2^10 

" + *4 _ir 2 \ 3 - V + 2 + 4 

L l( L J + L 2>10 , V 2 / 2 \ ( , * 2 loV L l - M 

4 + T" l”lO + V + + j 


- TT, 


’2 . L i( L i - 1 - L z) . 3 l 4 o . V: 


fh - l 2 \ 

" 1T 2 1T 3\ 4 y 


+ (rL 2 - L 1 )^“ 1T 4 L 2 ir 2^ ir 3 + U + \ 


L l +L 2 / * iO) +L /"10 


IT , 


+ TT , 


2 2 / 2\ 4 1 "2 


^ 2*10 


(1T3 + 1) + ^2^3 (^2 + 2) + 2 


“ TT, 


2 2 
L 2 L l + L 2 
4 + 4 + 8 


3 *10 * 2 , Al , 

+ 2 — (^3 + !) 


4 


(B13) 
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Normalizing factors . - Displacement magnitudes are normalized by 
division by the magnitude at zero rotational speed. This condition re- 
quires that 


ir 


4 


= 0 


(B15) 


Therefore, for support 1 



= 2[L (2r + 1) - L ] 

* 4=0 


and for support 2 


(B16) 


N„ 


1*4=° 


= 2[r(L 2 - L x ) - 2L 1 ] 


(B17) 


Setting each of the last two equations equal to zero gives the value of 
r for zero deflection at each support. From equation (B16) 


r N 1 =0 



(B18) 


From equation (B17) 



26 


r N 2 *0 



1 


Finally, equating these two equations to each other yields 



(B19) 


(B20) 


Physically, this set of conditions is impossible because although 
can be less than 1, 



Therefore, there is no set of conditions for which 


N 1 


N 2 

D 

I 

o 

1) 

D 


4 


=0 



TABLE I. - PARAMETER SELECTIONS 
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FIGURE 1. SKETCH OF SUPPORT SYSTEM SHOWING FORCE COMPONENTS 
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